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В. Г. МУРАДОВ, В. Я. ГАБЕСКИРИЯ 


О СОСТАВЕ ПАРА НАД ТВЕРДЫМИ СПЛАВАМИ, 
БЛИЗКИМИ ПО СОСТАВУ К СОЕДИНЕНИЯМ 
Са5! И Са>РЬ 


В работе одного из авторов [1] было сделано предпо- 
ложение, что при сублимации твердых бинарных сплавов 
Са, 51 и Са, РЬ, близких по составу к наиболее прочным 
соединениям Са5! и Са2РЬ, в газовую фазу переходят не 
только атомы компонентов сплава, но и молекулы этих 
соединений. 


Для проверки выдвинутой гипотезы нами исследован 
состав пара двух образцов Са, $1 (56 и 62 вес. % Са) и 
двух образцов Са, РЬ (66,5 и 70,8 ат. $ Са) в интервале 
температур от 400 до 850°С на масс-спектрометре МИ- 
1305. Испарение осуществлялось из цилиндрического тиг- 
ля, сделанного из нержавеющей стали 1Х18Н9Т, нагре- 
ваемого вольфрамовой печью сопротивления. Температу- 
ра определялась по току нагревателя печи, который гра- 
дуировался ‚специально на отдельном вакуумном посту 
с помощью хромель-алюмелевой термопары. Условия 
съемки масс-спектрограммы были следующие: напряже- 
ние ионизации менялось от 40 до 70 в, ток эмиссии като- 
да 1,5 ма, ускоряющее напряжение 2 кв. Во время опы- 
тов давление в области ионного источника было (2—3): 
. 10-6 тор, а в области анализатора (7—8) `10-7 тор. 

Расшифровка снятых масс-спектрограмм показала, 
что из твердых образцов Са, $1 в исследованном темпе- 
ратурном интервале сублимируют только атомы кальция. 
При этом отношение ионного тока атомов кальция к фо- 
ну составляло более 1000. Следовательно, пока можно 
утверждать, что если и присутствуют в газовой фазе мо- 
лекулы Са$1, то в количестве не более 0,1%. Однако сле- 
дует отметить, что проведенные эксперименты были до- 
статочно далеки от равновесных условий. Поэтому их 
пока можно рассматривать как ориентировочные. 

$ 


При изучении состава пара над твердыми образцами 
Са, РЬ уверенно удалось зафиксировать только ионные 


| (Са РБ) 
Г (РЬ) 


<0,01%. Отметим, что изученные образцы были частично 
окислены, так как на спектограмме были четко зафикси- 
рованы ионные пики, соответствующие окислам РЬО и 
РЬО.. Атомы кальция обнаружить не удалось, что кажет- 
ся очень странным, так как проведенный спектральный 
анализ конденсата показал в нем наличие большого ко- 
личества кальция. 

_ Таким образом. можно утверждать, что вопрос о суще- 
ствовании молекул. Са! и СазРЬ в газовой фазе остается 
пока открытым, так как четко выделить присущие им 
ионные пики на общем фоне в условиях наших опытов не 
удалось. 


Токи изотопов свинца. При этом отношение 
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ОБ ОЦЕНКЕ ТОЧНОСТИ ОТСЧЕТА ВРЕМЕНИ 
В ЭФФУЗИОННЫХ ОПЫТАХ БЕЗ ПРИМЕНЕНИЯ 
| ЗАТВОРА 


При измерении упругости паров по методу Кнудсена 
для точного отсчета времени эффузии обычно использу- 
ют различные конструкции затворов, перекрывающих мо- 
лекулярный пучок в процессе разогрева и охлаждения 
печи [1]. Такие устройства достаточно сложны в изготов- 
лении, хотя обеспечивают высокую точность отсчета вре- 
мени эффузии (порядка |—2% при экспозиции в 10 ми- 
нут). Покажем, что если использовать форсированный 
разогрев и быстрое охлаждение печи, обладающей не- 
большой тепловой инерцией, то можно обойтись без при- 
менения затворной системы. Для этого необходимо, что- 
бы систематическая ошибка, связанная с дополнитель- 
ным напылением исследуемого вещества на приемник в 
период разогреваьи охлаждения печи, была не больше 
общей ошибки эффузионного метода, составляющей 
обычно 15—20%. 

Рассматриваемая погрешность определяется несколь- 
кими независимыми факторами: 


|1) соотношением между временем разогрева и охлаж- 
дения печи (т!) и временем основной экспозиции (то), 
когда температура (То) поддерживается строго постоян- 
ной; 

2) скоростью разогрева и охлажления печи: 

3) природой вещества, т. е. зависимостью упругости 
паров исследуемого образца от температуры. 

Рассмотрим влияние каждого из этих факторов, ис- 
пользуя экспериментальные данные по упругости паров 
над системами Са—5$1, Са—РЬ и Са—5п [2]. При этом 
ограничимся рассмотрением только таких эксперимен- 
тальных точек, которые получены при малом времени эф- 
фузии то (10—30 мин.), так как легко понять, что с рос- 
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том то ошибка в отсче?те интервала времени эффузии бу- 
дет быстро уменьшаться. Интересующие нас эксперимен- 
тальные данные собраны в таблице 1. 


Таблица 4 


К расчету ошибки напыляемой на приемник массы 
в эффузионных опытах без применения затворной 
системы 
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Аз= 1 мии. 

ДлЯ Дн 

(11— 2,5) 

МИН. МИН. 
6,6 | 42 
3.8 1,9 
13.5 | 51 
147 | 63 
53 | 27 
104 | 559 
34 | 0,6 
7,2 4,6 
ии 
29 | 09 
2.9 0,9 
_ 59 й 1,8 
10.1 6.6 
4.8 | 1.1 
1.2 | 01 


В пятом столбце таблицы | приведены средние ско- 
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рости разогрева печи 





а }. измеренные в последние 
т 


две минуты перед достижением, заданной температуры То. 
Скорость охлаждения печи во всех опытах была порядка 
80—120 град/мин. 

Согласно формуле Кнудсена для эффузионного мето- 
да [1], систематическая ошибка, допускаемая при измере- 
нии количества напыленной массы на приемник без ис- 
пользования затвора, запишется следующим выраже- 
нием: 


ира 
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Шо Ро Со 


т: 


(0 


Даже если принять, что температура (Т) меняется в 
период разогрева и охлаждения печи со временем (т) по 
линейному закону, выполнить точное интегрирование в 
соотношении (1) представляется очень трудоемкой рабо- 
той. Однако, разбив время разогрева и охлаждения пе- 
чи т; на достаточно небольшие интервалы Ал, внутри 
которых упругость паров (Р) изменяется незначительно, 
можно приближенно рассчитать отношение (1), заменив 
точное интегрирование суммированием, т. е. 


Ш, Ве Ау Т. - Р; 
По “о Ро | == 1 ит, 


Анализ экспериментальных данных показал, что наи- 
более разумно выбрать Ал=| мин. Переход к ДАт= 


(2) 





т 
0,1 мин. уменьшает искомую величину у. на 20—954 
о 
(таблица 2), но резко увеличивает объем вычислитель- 
ных работ. 


Таблица 2 


Ш, 
Результаты расчета | —| при разных Ах для слу- 
| и 
чая № 12 из таблицы 1. 
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о /0 
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| 14,7 
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. т 
Результаты вычислений отношения Е по формуле 
т 


| о 
(2) при Алх= 1 мин. приведены в седьмом столбце табли- 


цы 1. В 8-м и 9-м столбцах этой же таблицы указаны 
погрешности, рассчитанные без учета 1,5 и 2,5 последних 
минут разогрева (охлаждения) печи перед достижением 
установившейся температуры То. 


„ [т 
Анализируя полученные результаты вычислений 
| та 


о 


для разных условий (таблица 1), можно прийти к сле- 
дующим выводам: 


п, \ оказывает 
1} Основное влияние на величину Е 


Шо 
т: | 
отношение Е. Значительно слабее исследуемая погреш- 
Со 
ность зависит от других факторов, т. е. скорости разогре- 
ва (охлаждения) печи и природы вещества. 

2) При выборе основного времени экспозиции (ло) не 
менее |5 мин. систематическая погрешность, вносимая в 
результаты эффузионных опытов, будет лежать в преде- 
лах случайной ошибки, т. е. не превышать 20%. 

3) Эта систематическая погрешность практическн 
полностью устраняется, если расчет упругости паров при 
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малых экспозициях (15—30 мин.) провести в предполо- 
жении, что эффузия при заданной температуре То про- 
исходила не то, а (10.-+2,5) мин. 

В заключение отметим, что исследованная ошибка бу- 
дет в действительности еще меньше за счет того, что в 
период разогрева печи упругость паров всегда меньше 
равновесной, принятой нами при расчетах. 
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ЭЛЕКТРОПРОВОДНОСТЬ БАРИЕВЫХ 
И СВИНЦОВОСИЛИКАТНЫХ СТЕКОЛ 


Большинство научных данных [1, 2, 3, 4] в настоящее 
‚время свидетельствуют о том, что при температурах ни- 
же 600°С силикатные стекла обладают преимуществен- 
но ионной проводимостью. Однако некоторые типы сте- 
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кол (бариевые, ванадиевые и некоторые другие) 65, 6; 
обнаруживают проводимость, свойственную электронным 
полупроводникам. Суждения о характере проводимости 
бесщелочных стекол до сих пор противоречивы. Так, на- 
пример, Штраус [7] утверждает ионную проводимость з 
свинцовосиликатных стеклах, обусловленную незначи- 
тельным количеством примеси щелочи. Павлова Г. А. 
[8] утверждает наличие электронной проводимости в 
этих стеклах. Незначительные (до 4 мол. %) примеси 
ионов щелочного металла в сложных стеклах [9] не мо- 
гут изменить электронный характер проводимости. 

В последнее время возрос ‘интерес к изучению сили- 
катных соединений: в кристаллической фазе.. Согласно 
ряду авторов [10,11], соединения в аморфной и кристал- 
лической фазах имеют ряд отличительных особенностей: 
величина проводимости в кристаллической фазе для ще- 
лочных и некоторых других силикатов значительно ни- 
же, чем в аморфной, а влияние состава на электропро- 
водность слабее выражено в кристаллической фазе. Кро- 
ме того, было установлено [10], что характер кристал- 
лизации (величина размеров кристаллов) мало сказыва- 
ется на-величине электропроводности. Все это говорит в 
пользу большей роли электронной проводимости в кри- 
сталлическом веществе. `- 

В настоящей работе исследовалась зависимость объ- 
емной проводимости от температуры в бариевых и свин- 
цовосиликатных стеклах. Бариевые соединения (28 мол. 
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% ВаО) были предварительно закристаллизованы, свин- 
цовые стекла исследовались в аморфной фазе. Образцы, 
предназначенные для исследований, подвергались испы- 
таниям на «формовку». Было установлено, что до темпе- 
ратуры 450°С и плотности тока —10-5 а/см? сопротивле- 
ние их в постоянном режиме не зависит от времени. Вы- 
полняемость закона Ома и независимость сопротивления 
от времени дополнительно подтверждают преобладание 
электронной проводимости указанных стекол. Измерения 
велись по обычной схеме с охранным кольцом компенса- 
ционным методом. В качестве индикатора в импульсном 
режиме использовался осциллограф ИНО-1 (С1-—4). В 
постоянном режиме использовался зеркальный гальвано- 
метр. Импульсы прямоугольной формы имели частоту 
следования 25 гц и напряжение | кв. 

Исследуя зависимость проводимости ВаО`$1О0. от тем- 
пературы в импульсном режиме (рис. 1), легко видеть, 
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что характер этой зависимости различен для различных 
времен. Меньше всего температура сказывается на ма- 
лых временах (2`10-”сек.), что указывает на слабую за- 
висимость быстрых состояний от температуры. 
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Рис. 2 


Интересна зависимость |2 1,= Р (1/Т) для свинцового 
стекла (рис. 2). Легко видеть, что график распадается на 
2 прямолинейных участка, соответствующих энергиям ак- 
тивации 0,65 эв и 0,69 эв. Можно предполагать, что в дан- 
ном случае освобождались лишь неглубоко расположен- 
ные состояния. Не исключена возможность, что при более 
высоких температурах обнаружатся и более глубокие 
энергетические: состояния. 
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_В. А. ГУБИЧЕВ 


ВОЛНОВЫЕ СВОЙСТВА ЧАСТИЦ И УРАВНЕНИЕ 
ШРЕДИНГЕРА В КУРСЕ ОБЩЕЙ ФИЗИКИ` 


Волновые свойства частиц так же, как корпускуляр- 
ные свойства волн (полей), представляют наиболее фун- 
даментальные свойства материи. Не во всех вопросах 
этой проблемы в современной физике достигнута доста- 
точно предельная ясность, но основы проблемы разрабо- 
таны достаточно хорошо и нашей задачей в курсе общей 
физики является выработка отчетливых представлений 
об этих основах в сознании обучающихся. 

Основным законом, описывающим волновые свойст- 
ства частиц, является уравнение Шредингера. Это опыт- 
ный закон, играющий в данном разделе физики такую 
же роль, какую второй закон Ньютона играет в динами- 
ке. Доказательством его справедливости является совпа; 
дение выведенных из него следствий с опытными фак- 
тами. | 

Собственно, с этого утверждения и начинаются труд- 
ности. при изучении уравнения Шредингера. Вся подгото- 
вительная работа, заключающаяся в описании догадки 
Луи-д6-Бройля о волновых свойствах частиц и опытов, 
подтверждающих эту догадку, воспринимается с довери- 
ем, но когда записывается уравнение Шредингера и объ- 
является, что это опытное уравнение, то доверию прихо- 
дит конец и никакая апелляция ко второму закону Нью- 
тона здесь не помогает. Действительно, нетрудно пред- 
ставить себе опытное происхождение такого простого со- 
отношения, как второй закон Динамики, и совершенно 
невозможно представить это по отношению к уравнению 
Шредингера — этому сложному уравнению. 

Положение усугубляется еще и тем, что, как правило, 
студенты к этому времени не ознакомлены с волновым 
уравнением в курсе математики. Не способствует отчет- 
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ливому усвоению идей и то, что далеко не во всех учеб- 
никах и, соответственно, лекциях четко разграничено из- 
вестное и неизвестное о% — функции — этой загадочной 
функции, которая характеризует состояние частицы, а 
сама остается неведомой. У студента при этом часто 
остается впечатление, что только ему по какой-то при- 
чине остается неясным смысл этой функции, ‘а другим это 
известно. Появляется неверие в свои силы. 

Зная об этих трудностях, некоторые преподаватели 
пединститутов, даже на физических отделениях, при изу- 
чении корпускулярно-волнового дуализма материи огра- 
ничиваются только описанием опытов, тем более, что зна- 
ние ф — функции и уравнения Шредингера не предусмат- 
ривается программой. Удивительное дело! Студенты тех- 
нических вузов должны знать об этом, а вот в програм- 
ме курса общей физики на физических отделениях пед- 
институтов знакомство с раВневив Шредингера не 
предусмотрено! 

Не подлежит, однако, сомнению, что общие принципы 
корпускулярно-волновой тебрии столь важны, что пред- 
ставляется совершенно необходимым более широкое зна- 
комство с элементами этой теории еще в курсе общей фи- 
зики. А указанных выше трудностей можно избежать, 
если не записывать сразу уравнение Шредингера в об- 
щем виде, заявляя, что это волновое уравнение, а пока- 
зать, что оно является волновым, опираясь на знания о 
волнах, полученных в предыдущих разделах курса физи- 
ки. При этом рекомендуется придерживаться следующего 
плана: 

а) Описание гипотезы де-Бройля о волновых свойст- 


| и р 
вах частиц и обоснование формулы де-Бройля ).=-= о) 


на основе пропорциональности массы и энергии. 

6) Эксперимёнтальное подтверждение — волновых 
свойств частиц в опытах Девиссона и Джермера, Тарта- 
ковского и др. Не забыть сказать, что из этих опытов 
можно определить длину волны частицы и что она сов- 
падает с длиной волны, предсказанной де-Бройлем, что 
является еще одним подтверждением справедливости 
его идей. 

в) (Основной пункт, из-за него и статья писалась). 
Если в механике, электромагнетизме при исследовании 
волн известно, какая величина совершает колебания: в 


г 


= 
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звуковой волне — это давление воздуха, в электромаг- 
нитной — напряженности электрического и магнитного 
нолей и т. д., то природа волнового процесса, сопровож- 
дающего движение частиц, неизвестна. Какая физиче- 
ская величина здесь совершает колебания неизвестно. 
Однако опыты показывают наличие дифракционных мак- 
симумов и минимумов, значит есть какая-то волна и 
есть величина, совершающая колебания. Эту неизвестную 
величину обозначили буквой + и назвали волновой фун- 
кцией или + — функцией. 

Все волны, распространяющиеся в одном направле- 
нии, описываются одинаковыми уравнениями. Для ме- 
ханической волны, распространяющейся вдоль оси х, та- 
ким уравнением будет: 

х 
у- ут [ =). 
у 
для электромагнитной: м ' 


5 


Х Хх 
В, иь (1%). НЕН,Зть (1 }. 
у А 


По аналогии, волновое уравнение для частицы, движу- 
щейся вдоль оси х, запишется в виде: 


! <] | г 
ф —= Фо © о | 
Дважды дифференцируя это уравнение по х, получаем: 


д:ф 0 х 
= уз шо = 





или 
2 2 
2 шв 


а ти 


Это волновое уравнение для движущейся частицы, 
т. е. уравнение Шредингера, его легко привести к обычно- 
му виду, если учесть, что: 











(0 Эх Ю Е 
А = а тогда 
22 4к _ 412 32 у? Вт? т | 
ее - Ех. 


Ек — кинетическая энергия частицы. 
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После подстановки последней формулы уравнение 
волны приобретает вид: 
д 8 = т о 
АРА 
‚9х? т. 18 


Это уравнение Шредингера для стационарного одно- 
мерного движения свободной частицы. 

Если частица движется в силовом поле, то ее общая 
энергия Е есть сумма потенциальной и кинетической 
энергии и следовательно: 


а 


где (О — нотенциальная энергия. 

Произведя эту подстановку, получаем уравнение Шре- 
дингера в более общем случае: 

0?Ф` Вт: т 
Зе 

После получения этого уравнения совершенно естест- 
венно воспринимается обобщение его для случая части- 
цы, движущейся произвольно: 

92ф оф 0% Вет | 
дх? т ду? у 02 Г 1? о м 

Приведенные рассуждения не следует рассматривать 
как вывод уравнения — это преобразование волнового 
уравнения, данного в интегральном виде (данного опыт- 
ным путем), к дифференциальному виду. Такой подход 
облегчает недостаточно эрудированному в математике 
студенту понять смысл уравнения Шредингера, а в даль- 
нейшем помогает и в выявлении физического смысла ф — 
функции. | 

д) Далее традиционным методом (поэтому он не 
описывается подробно), путем сравнения оптических 
дифракционных картин и картин дифракции электронов 
можно показать, что, хотя физический смысл самой %— 
функции неизвестен, квадрат модуля Ф — функции пред- 
ставляет вероятность нахождения частицы в месте с ко- 
ординатами х, у, 7 в момент времени {, показывая таким 
образом, что Ф — функция имеет статистический харак- 
тер. 
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Н. П. ГРИФЦОВ 


НЕКОТОРЫЕ ВОПРОСЫ ДИНАМИКИ ТОЧКИ 
ПЕРЕМЕННОЙ МАССЫ В ЭЛЕМЕНТАРНОМ 
ИЗЛОЖЕНИИ 


В наш век — век основания космоса большой интерес 
учащихся вызывает динамика полета ракеты. Однако 
предусмотренный прбграммой средней школы материал 
не дает возможности провести соответствующие, хотя бы 
приближенные расчеты, так как они связаны со специ- 
альным разделом механики — динамикой точки пере- 
менной массы, не проходимым в школе. 

Этот раздел довольно труден, но отдельные важные 
вопросы нам представляется возможным рассматривать 
во внеклассной работе, опираясь на знания учащихся по 
физике и математике, не выходящие за рамки школьнон 
программы. 

Данная статья представляет собой несколько перера- 
ботанную и сокращенную лекцию, прочитанную нами в 
школе юных физиков при Ульяновском пединституте. 


1. Понятие точки переменной массы 


В механике Ньютона массы движущихся тел считают- 
ся при их движении постоянными (в некоторых случаях 
рассматривается ‘изменение массы скачком на некоторую 
величину. Примером может служить выстрел из орудия, 
когда масса последнего уменьшается на величину массы 
снаряда). 

Однако в природе, технике многочисленны случаи, 
когда масса движущихся тел меняется непрерывно вслед- 
ствие отделения от них частиц, или, наоборот, присоеди- 
нения. 

Такие тела называют телами переменной массы. Для 
них масса является непрерывной функцией времени: 


т=ш (1) 


Приведем несколько примеров таких тел 


1. При работе магнитофона масса‘одной кассеты не- 
прерывно уменьшается, а другой — увеличивается. 

2. Масса айсберга, плавающего в море, уменьшается 
по мере таяния льда или увеличивается при его намора- 
живании. 

3. При поднятии якоря масса свисающего участка 
цепи непрерывно уменьшается, при опускании — увели- 
чивается. | 

4. Масса ракеты при работе двигателей непрерывно 
уменьшается за счет охбрасывания продуктов- сгорания. 

Можно привести целый ряд примеров движения тел 
переменной массы. Даже наша Земля, вследствие паде- 
ния на нее метеоритов, непрерывно увеличивает свою 
массу. | 
Впервые уравнения движения тел переменной массы 
были получены выдающимся русским ученым И. В. Ме- 
щерским в 1897 г. Им фактически был создан новый раз- 
дел механики — теория движения тел переменной мас- 
сы, являющейся научной основой современной ракетоди- 
намики.. 

Мы будем рассматривать точку переменной массы, 
под которой будем понимать материальную точку, в ко- 
торой сосредоточена масса всего тела, изменяющаяся не- 
прерывно по некоторому закону е 


ш=иш ({) 


Любое тело можно рассматривать как материальную 
точку, если его размерами по сравнению с прохолимыми 
расстояниями можно пренебречь. 


2. Уравнение движения точки переменной массы 
(уравнение Мещерского) 


Перед выводом уравнения необходимо ознакомить 
учащихся с теоремой об изменении количества движения 
механической системы, которая имеется в любом учебни- 
ке физики. Эту теорему можно сформулировать следую- 
щим образом: изменение количества движения системы 
за некоторый промежуток времени равно сумме импуль- 
сов всех действующих на тела системы внешних сил за 
тот же промежуток времени. 
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ЕО 
Математическая запись: 
Если обозначим количество движения системы в на- 
чальный момент ‘времени (Фо, а через время А1{—0Оу, то 


,—9,=ХРА! (1) 
Применим эту теорему к выводу уравнения Мещерского. 


Пусть точка, масса которой изменяется с течением 
времени, в момент времени %=0 имела массу т и ско- 


рость У по отношению к Земле. Во время движения на 


нее действует результирующая внешних сил Е’\". 
За время 1 к этой точке присоединяется частица 
массы Ат, имеющая в момент присоединения скорость 


по отношению к Земле (Ц и на нее действует результи- 


рующая внешних сил Е.”. 

Будем рассматривать эти две массы как одну систе- 
му. Тогда масса системы не меняется, и можем приме- 
нить теорему (1): 


Количество движения этих масс в момент %=0: 
о = шУ -- А ШП 
После присоединения массы Ат скорость движущей- 


ся точки изменится на величину ЛУ и обе массы будут 
двигаться уже с общей скоростью У-+ АУ. 


Количество движения в этом случае: 
О, = (ш+ Дм) (У АУ) 


Применяя теорему (Г), получим: 


У — ЛшУ + шАУ -- АМАУ — шУ — АшИ == 
= (РЕ) ЛЬ. 
учитывая, что присоединяемая масса Аш, а также изме- 


нение скорости АУ величины малые, можем пренебречь 


их произведением. Приведя подобные члены шУ и —ш\У 
и сгруппировав: слагаемые, содержащие Ат, получим, 
обозначив сумму внешних сил одной буквой Е‘. 


(У— 0) Дм + мАУ = 2 д! 
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Поделим все члены полученного уравнения на А{и 
перенесем слагаемое, содержащее Ат в правую часть. 
В результате получим искомое уравнение Мещерского: 

^ ее 2 нба ИИ а 
п к 2 
м р Их (2) 

Однако, обычно его записывают в несколько другом, 
более удобном для применения виде, вводя так называе- 
мую относительную ^ скорость Ч, представляющую со- 
бой скорость присоединяемой массы /ААт относительно 
массы т. 

Очевидно, что 9—У=:(.. 

После введения относительной скорости Ц. получа- 
ем окончательно: 


Жи Ат 
д № (3) 


Более строго это уравнение получится путем перехода к 
пределу при АЕ -— 0, т. е. 
Е Е Е 
р, шт —=—— 
В этом случае получим уравнение Мещерского в 
обычно рассматриваемом виде 
ЧУ Ен т 
ЕЕ ы 4 
Однако для применения к расчету реактивной силы 
достаточно и полученной нами элементарной формы. 





= Р°-- Оо 





3. Анализ уравнения Мещерского 


1. Выясним размерность второго слагаемого правой 
части: в технической системе единиц получаем его раз- 
мерность в кГ, в системе СИ — в ньютонах. Таким обра- 


.  АШ 
зом член (Ц А имеет размерность силы и называ- 
ется реактивной силой К, т. е. 
Вы (5) 
= 5 
о 1 
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2. Рассмотрим, как направлено ускорение, создавае- 
мое реактивной силой в случае присоединения или отбра- 
сывания массы. Для удобства анализа примем, что внещ- 
ние силы отсутствуют. В этом случае уравнение Мещер- 
скогс можно записать в виде 


и Ат 
ша = Ч, 
Ау а — 
так как -„__ представляет собой ускорение а. 
а) К движущейся точке массы т присоединяется ча- 
Ат _ Ат 








стица Аш. В этом случае МОЖНО 


д: 
рассматривать как скорость изменения массы. 


Тогда 0, иа совпадают ‘по направлению. 
6) От движущейся точки отделяется масса Ат. Сле- 


Ат 
довательно, < ) и ускорение направлено в сторону, 


обратную скорости движения частиц относительно мас- 
сы м. 

Этот случай наиболее интересен, так как имеет место 
при движении ракеты. В дальнейшем мы и будем его. 
рассматривать. 

3. Если масса движущегося тела не меняется, то 


Ат 
АЕ. 
ный 2-Й закон Ньютона: 

МУ = 
А —) . 


—= О и уравнение Мещерского переходит в обыч- 





Ш 


4. Применение уравнения Мещерского к Орояснению 
полета ракеты 


Запишем уравнение (3) в виде 


А. 
ко е 6! 
. Е =Р Е. р (61, где 
К — реактивная сила, определяемая формулой 
те г т 
К = % —- Из уравнения (6) следует, что даже 
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если ракета движется вне поля внешних сил, т. е. Е° == 0, 
она получает ускорение а за счет реактивной силы, зави- 
сящей только от относительной скорости истечения про- 
дуктов сгорания и скорости изменения массы. 

Для современных ракет Ч. имеет порядок 2—3 
км/сек (для химических топлив) и при расчетах может 
приниматься постоянной (для установившегося режима 
работы двигателей). | 





т 
Скорость уменьшения массы, т. е. 2. при этих ус- 
ловиях также можно считать постоянной и ие обыч- 
х\ 


но называют секундным расходом топлива М. Число М 
показывает, какое количество топлива сгорает при работе 
двигателей в 1 секунду. 

Таким образом, величина реактивной. силы целиком 
определяется скоростью истечения продуктов сгорания 


НЕ ы Л 
Ч и секундным расходом топлива М -= г 


Рассмотрим характерный пример, из которого видно, 
какое огромное количество топлива сжигается современ- 
ными ракетами для получения необходимой реактивной 
СИЛЫ. | | 

В журнале «Авиация и космонавтика» № Зза 1965 г. 
опубликованы некоторые данные об американских стра- 
тегических ракетах. Одна из них, межконтинентальная 
баллистическая ракета «Титан», имеет длину 35 м, диа- 
метр до 3,5 м. Вес ее на старте 100 тонн, суммарная тяга 
двигателей 160 тонн (реактивная сила). 

Рассчитаем секундный расход топлива этой ракеты, 


приняв относительную скорость истечения газов Шо ==2,5 
км /сек. 
Дано: 


Ю = 160000 кГ =. 1570000 н. 


Ч = 2500 м/сек. 





Аш 


Определить М == 5. 
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Решение: 

















= сы 
Применяем формулу реактивной силы В = (о ее | 
6 Аш _ К 
тсюда их 0. 
Подставляя данные, получим: 
р кг м 
лш 1570000 ИЕ Е 
> 2900 
сек 


Для сравнения можем привести некоторые данные наше- 
го космического корабля «Восход» (данные были опубли- 
кованы в газетах). 

Корабль «Восход» имел 7 двигателей суммарной тя- 
гой 650 тонн. Вес корабля без учета веса ракет-носителей 
5320 кГ (о стартовом весе данных в печати не появля- 
лось). | 

Однако при расчетах следует иметь в виду, что не все 
двигатели работают одновременно, так как для получе- 
ния необходимых скоростей применяют многоступенча- 
тые ракеты, впервые предложенные нашим замечатель- 
ным ученым К; Э. Циолковским. 
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ОРГАНИЗАЦИЯ И ПРОВЕДЕНИЕ 
СПЕЦПРАКТИКУМА ПО ФИЗИКЕ 
В ПЕДИНСТИТУТЕ 


Спецпрактикум по физике для студентов старших 
курсов физико-математического факультета отделения 
физики пединститута является одной из завершающих 
практических дисциплин, способствующих подготовке 
учителей-физиков к практическому осуществлению задач 
политехнического обучения учащихся средней школы: 
Проведение этого практикума, по которому отсутствуют 
программы, должно базироваться на достижениях совре- 
менной физики и техники, современных методов физи- 
ческого исследования, и соответствовать той роли, какую 
имеет физический эксперимент в преподавании физики в 
средней школе и физические измерения, являющиеся ос- 
новой для понимания количественных законов физики. 


Нами был поставлен спецпрактикум по полупроводни- 
кам и полупроволниковым приборам, которые в настоя- 
щее время занимают центральное место в физике твердо- 
го тела, являясь основой технического. прогресса во мно- 
гих отраслях техники. Организация такого спецпракти- 
кума позволяет углубить и дополнить теоретические све- 
дения, получаемые студентами в лекционном спецкурсе 
по полупроводникам, изучить методы исследования и из- 
мерения физических параметров полупроволниковых ма- 
териалов и полупроводниковых приборов, изучить прин- 
ципы создания установок для проведения эксперимен- 
тальных исследований. „ 


Нами поставлены следующие лабораторные работы: 


1. Изучение полупроводниковых термочувствительных 
сопротивлений. ь 
2. Изучение термоэлектрических явлений. 
`3. Изучение эффекта Холла в полупроводниках. 


4. Изучение выпрямительных свойств гермениевых и 
кремниевых диодов, стабилитронов, туннельных диодов. 

5. Изучение фотоэлектрических свойств полупровсл- 
НИКОВ. 

6. Измерение удельного сопротивления полупровод- 
ника компенсационным методом и изучение распреде- 
ления потенциала вдоль полупроводника. 

7. Изучение работы полупроводникового триола. 

8. Изучение полупроводников в школе. 

9. Изучение сегнетоэлектриков. 

10. Изучение ферритов. 

11. Изучение работы фотодиодов и фототриодов. 

12. Определение контактной разности потенциалов. 

13. Изучение нелинейных полупроводниковых сопро- 
тивлений (варисторов). 

В руководствах к лабораторным работам указывается 
литература для изучения теории и методики работы, 
программа работы, описание экспериментальной установ- 
ки, методика проведения эксперимента, контрольные во- 
просы. К каждой лабораторной работе составлены прог- 
раммы для обучающей машины «Ласточка». 

Проведение спецпрактикума по физике имеет ряд сво- 
их особенностей по сравнению с практикумом по общей 
физике: | 

1. Теория вопроса подробно изучается самостоятель- 
но по указанной литературе. Студент должен научиться 
самостоятельно из рекомендованной литературы выбрать 
необходимый ему материал по данной теме работы. 

2. В руководствах к работе не приводятся таблицы 
измерений и формулы ошибок. Студент должен самосто- 
ятельно составить таблицу измерений и оценить точность 
и реальность полученного экспериментального резуль- 
тата. 


3. В руководствах предлагается произвести некоторые 
измерения разными методами, сравнить полученные ре- 
зультаты, выбрать наиболее рациональный и более точ- 
ный метод измерения физических величин. , 

4. Большинство измерительных технических приборов 
и установок студенты самостоятельно осваивают соглас- 
но’заводским инструкциям по их эксплуатации. 

Лабораторные работы, как правило. выполняются ин- 
дивидуально по графику цикличной системы, т. е. студент 
выполняет работы в круговом порядке. Занятия прово- 
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дятся по расписанию деканата или по свободному распи- 
санию, когда каждхи студент выбираст ему удобный 
день занятий в лаборатории и посещает лабораторию в 
этот день. 

В лаборатории имеется почти вся справочная литера- 
тура, научная и методическая литература по тематике 
лабораторных габот. 

Студенты готовятся к выполнению лабораторных ра- 
бот как дома, так и в.лаборатории до выполнения работы. 

При выполнении работы студент оформляет отчет 
(название работы, основные рабочие и принципиальные 
схемы, полученные результаты измерений и вычислении с 
оценкой точности полученных измерёний и вычислений). 

При приеме отчета по работе обращается . внимание 
на следующее: | 

1. Знание физического смысла измеряемой величины, 
методики ее измерения, теории, на которой базируется 
работа. 

2. Физические ‘принципы устройства измерительного 
прибора, его принципиальная схема, умение пользо- 
ваться им. 

3. Знание степени точности результата измерений и 
вычислений. 

4. Применение изученного явления в науке, ‘технике, 
народном хозяйстве. 

В лаборатории спецпрактикума по полупроводникам 
используется обучающая машина «Ласточка». К ней со- 
ставлены программы для всех работ. 


`Машина используется в основном в режиме «кон- 
троль» для проверки подготовленности к выполнению ла- 
бораторных работ, приема зачетов. Для слабо подготов- 
ленных студентов машина иногда включается в режим 
«самоконтроль», т. е. как обучающая. 


Программы к каждой работе включают материал, 
который студенты должны самостоятельно проработать 
по рекомендованной литературе. Он располагается в ло- 
гической последовательности и делится на небольшие 
порции. Контроль за усвоением каждой порции и осуще- 
ствлялся нами на обучающей машине. 

Применение таких программ для проверки усвоения 
студентами рекомендованной литературы стимулирует их 
самостоятельную работу, выявляет глубину изученного 
материала, сокращает время преподавателя. Преподава- 


— 


тель может больше времени уделить обучению студентов 
правильной методике проведения эксперимента и получе- 
нию наиболее точных результатов, более глубокой про- 
верке знаний студентов. 

Однако применение обучающей машины ни в коем 
случае не исключает собеседования по отдельным вопро- 
сам методики выполнения эксперимента, обоснованию 
ответов, полученных результатов, обобщению изученного 
материала. 


Э. М. БАШИРОВ` 


ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ПРЯМОУРОЛЬНОЙ 
ПОТЕНЦИАЛЬНОЙ ЯМЫ 


При расчетах соударений атомов часто [2] для просто- 
ты потенциал задается в виде 


> (<) 
У (г) = | —е (а<г-< \а) 
0 (г>> ›а) 


Цель настоящей статьи: определить параметры &, а и 
у по экспериментальным данным рассеяния, где & — глу- 
бина ямы, а — радиус отталкивающей стенки и у — чис- 
ло, показывающее во сколько раз радиус области притя- 
жения больше радиуса области отталкивания. 

Предположим, что задачу можно решать классически. 


9> 34 
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На рисунке схематично показано, как происходит рас- 
сеяние. При о >> уа частица не рассеивается (классичес. 
кая траектория представляет собой прямую). При < 
о < 14 имеет место рассеяние в диапазоне углов 
(о, 8,). При дальнейшем уменьшении прицельного пзё- 
раметра р частица натыкается на отталкивающий барь- 
ери не может рассеяться на угол больший, нежели $, — 
так называемый радужный угол. |1] Траектория частицы 
разгибается. Проходится диапазон углов (о, $, ). Рассея- 
ние на о’ происходит при некотором р-==рь. При изме- 
нении прицельного параметра в промежуткео < р < 
частица проходит диапазон углов (ол). Таким образом 
область углов от 0 до \, проходится три раза, аот 3, 
дол — один раз. 


Угол 2Ф, который описывает следящий вектор, когда 
частица движется из бесконечности и вновь удаляется в 
бесконечность, для различных диапазонов из-за нали- 
чия двух областей задания потенциала записывается по- 
разному, а именно: 





оС УД, 
? рАг я оЧг 
ф/ — Е ити СС ЗА ря 5 — 
Е ы 2 - р 
т и 1 — 2 Г и Е 2 
УД Го 








$э = == Бо т За 
и 2 и 
Е 


< 
> 
® 


Го — корень подкоренного выражения в [2; он равен 


1 


= р 
Г == (1 — Е) 2 . Отсюда можно найти ир, при усло- 
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1 
ВНИ, ЧТО Го = 1, то есть р; == 4 (1 - = 4 . Интегралы вы- 


числяются элементарно. Для Ф имеем: 











: о в Ра ь 

®, — агс$1а _- -- агссоз в (5; < в.<»а), (1). 
УС, УС. 

т ри * 

$> == агсзт`,, -- агссо$ о ЕВЕ: — в (0- < в,), 
(2) 
Е | 

где в Е и Е — энергия налетающей частицы. От- 


сюда легко получить выражение для радужного угла 
рассеяния %,, подставив в (1) (или в (2) значение р; и 


5 
о ’ 
| 1 


Ире -|- 2 170005 == агсс08 = 
Уи и У 





учитывая, что в области притяжения $ == 





(3) 


Из (3) ясно видна зависимость радужного ‘угла от энер- 
гии: с увеличением энергии налетающей частицы 3, 
уменьшается. 

Дифференциальное эффективное сечение для любых 9 
трудно выразить в аналитическом виде (из-за сложной 


зависимости р от % для О < рр). Однако удается 
его записать для очень больших углов рассеяния. 


Перенося в (2) первый член правой части. -налево и 
беря с0$ от обеих частей получившегося выражения, по- 
сле несложных преобразований получаем следующее 


уравнение относительно х — я 


УД, 
Черлхе — 4 (вв Ох нь м НТ 2: 
- (У? -- У + 1) с0з 2$] х 
-- 4урзшф (1 -- созф)х? — 2 102 (У? 6 ь + 1) х2 + 
-- $110 ф = О (4) 
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Рассмотрим область малых р, т. е. рассеяние близкое к 
180°. В этой области х << 1. Пренебрежем членами с 
х" |п>2| по сравнению со свободным членом. Анализ 
показывает, что это приближение проходит для углов 
рассеяния от 1177° до 180°. | 


| 
(Анализ проведен для у—9, в — ыы 
Таким образом для углов, рассеяния от 177” до 180°, что 


соответствует прицельным' параметрам, лежащим в обла- 
сти О%р< 0.02 а, имеет место уравнение 


2х2 (Уи -Ны- 1) = $112, 





9 
откуда (учитывая, что для отталкивания ф = —_) 
с Ч? с У? (12 31п $ на 








(5) 
Добавляя к уравнениям (3) и (5) еще выражение для 
полного сечения, которое для прямоугольной потенциаль- 
ной ямы будет, очевидно, равно | 
о== (0), (6) 
получим систему трех уравнений (3), (5) и (6), решая 
которую относительно = ‚, би у, определяем искомые па- 
раметры. 
Из-за важной роли радужного угла в предлагаемой 
методике энергия налетающей частицы должна быть до- 
статочно мала [Е < О: 1е\|. - 


В заключение выражаю глубокую благодарность 
Л. Г. Яковлёву за помощь в работе. 
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СВЯЗЬ ПАРАМЕТРОВ СВЯЗАННОЙ РЕАКТОРНОЙ 
СИСТЕМЫ С ДЕТАЛЬНЫМ РАСПРЕДЕЛЕНИЕМ 
НЕЙТРОННОГО ПОТОКА` 


Запишем двухгрупповые уравнения стационарного со- 
стояния отдельной сборки в виде: 


м 5 \/ 9; — Хо $5 -- Хю $; == 0, (1) 
\/ У Др М $: — Хо ФЕ Кио Жо Фз = 0. 


Здесь Б, У и к—функции координат, значок 0 означает, 
что эти параметры не возмущены. Значок $ относится к 
величинам, связанным с описанием тепловой группы 
нейтронов, Г — группы быстрых нейтронов. Смысл симво- 
лов О, %, киф очевиден (1). 

Под сборкой понимается ансамбль из размножающе- 
го нейтроны материала. 

Пусть в сборку внесено возмущение путем малых из- 
менений одного или нескольких параметров и одновре- 
менного сближения данной сборки с некоторой другой 
критической сборкой. 

Изменения параметров. могут иметь разные знаки и 
должны быть такими, чтобы система никогда не прибли- 
жалась к состоянию мгновенной критичности. Сборка с 
возмущенными параметрами будет описываться уравне- 
НИЯМИ: 


1 0 
УЛ, Уф, — Хр, + ХЕ фе- $ (г) = р = (3) 


УГ: у: — ХЕ (1 - ВК, 93 + 





сне бе к (1) 
1 У; Я ' 
дс. 
^^, С; - к, №, В ф: = — ® (5) 
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Здесь фз; — 9 (Г, Ё), Ч = ф: (г, ), и. а 


В. = О, {г} + 80, (г*, Обр = 0 (г) -- 80, С), 
р ВЕ ма (г) но рае р —: и в (г) + бу, (г), (6) 
К. м - Ко] ое ТГ 0 (к: ра 


Взаимодействие, в общем случае, осуществляется как 
на тепловых, так и на быстрых нейтронах. Соотношение 
между числом тепловых и быстрых нейтронов, приходя- 
щих от второй сборки, зависит от характера среды, раз- 
деляющей сборки, и расстояния между сборками. 

Поскольку экспериментально наблюдается единый 
установившийся период для всей системы, то поступаем 
так, как будто переменные разделяются: 


: (Г, @) = $’ (г)Е (0), 
— Хх / 4 
ф (Г, 1) т (г) Е, (1), (7) 
С! (г, 1) = С: (г) Е (©. 

Считаем также, что вторая` сборка создает в первой 
нейтронные источники, распределение которых описыва- 
ется функцией -5 .(г) для тепловых нейтронов и функцией 
$ Кг) для быстрых, причем во времени интенсивность ис- 


точников меняется так же, как и мощность во второй 
сборке, но с некоторым запаздыванием т ‚ т. е. 


Ре (8) 
$: (г,) = $: (6) Е. (1— *,. (9) 


Запаздывание связано с распространением нейтрон- 
ной волны от одной сборки к другой и различно для бы- 
строй и тепловой групп нейтронов [2]. 





Нам понадобятся также уравнения для сопряженных 
потоков ФФ"; и Ф": 
УР 9", —— Хоф" ко ХУ ф "= 0, (10) 
УП У $: * — Хюфг" -- Хю$." = 0. (11) 
а) Умножим теперь уравнение (3) на $.*, а (10) на 


$. Г:(@), проинтегрируем по объему первой сборки и 
вычтем друг из друга. 


3.1 


ИХ 


6) Умножим также (4) на ф/“,а (11) на $, Р, (1, 
проинтегрируем по объему первой сборки и вычтем друг 
из друга. 


Сложим теперь результаты операций (а) и (б) и раз- 
делим на 


1 кг Уз фе? $. ЧУ. 
У! 


В результате получим: 











ав, (9 м8 И во 
о Е 
91 И | Г 
| (12) 
ав (1 
т == — № (0 - . мт (13) 


= 18, Уз) фе — 8 Уф. — ЗУ фг -- ЗУ ег — 
У, 





(14) 
. ^ и ИЕ / 25 ее. 3: 
— 0. УФ У Ф.. — 9): \/ Фг- \/ $ ] Чу | | К 5. фг" $: У, 
а 
СУ а. — У: - № 
(15) 
В = В, ] КУР фь ЧУ | к Хз фе" $5 ЧУ, (16) 
У! | 
— в 
= и где (17) 
в, (г) 
с = ©01п${; 
К, Уз Фь (Г) 
из = [[$. (г) фу Чу [$ : (г) фе ау] 
У! У! 
[к Уф" 9, Чу, (18) 
У! 


из = [3 [$ (Г) фз* ау -- м [$1 (г) фг* 9\]/ 
У! У! ` 


| | $5 (г) 95* Чу -- [$ (г) 9 ыы (19) 


У, У, 


Конкретная форма функций $. (г) и $;(г) зависит от 
выбора ядер замедления и диффузии. 

Повторяя аналогичные выкладки для второй сборки, 
можно получить вторую пару уравнений. Вводя привыч- 
ные обозначения п, (1) > Е, (®, С (> (8, Н->1*; 
3, — 8,*", запишем систему уравнений кинетики двойной 


связанной системы в окончательном виде: 


—9 ь 
о АН п; (#1 -- ХА Сы (0 +: И (4 — <>), 
{0 1. | 
ЧП. 4 2) = Но | > 
о и п. (В > Сы (| а п; (Е — 21), 





(20) 
ЧС; ). : 
Е АОИ ЗЕ “- п, (1), 
аСь (# . ь 
БЕ р ЕЕ 
1 Е Ср 
Здесь п; (1) и п. (1) — средняя плотность нейтронов 


(или общая мощность сборок) в первой и второй сборках 
соответственно. 

Си ({) и С. (® — средние плотности (или средние 
мощности источников) предшественников запаздывающих 
нейтронов 1-й группы в первой и второй сборках, соответ- 
ственно. 

12 И >, — коэффициенты связи; в общем случае 
112 5= 121. 

Ни р — среднее время цикла м. МГНО- 
венных нейтронов. 


В, и В, — эффективная доля запаздывающих 
а. 
Заметим, что знаменатель в (14) — (18) может быть 


выбран произвольно. Используемая здесь форма знаме- 
нателя выбрана таким образом, что при 1 —› О получается 
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обычное определение реактивности, например, (14) дает 
о — бк,/к, [3]. , | 
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Е. Д. ВОРОБЬЕВ, Е. Г. КАЛАШНИКОВ 


РЕШЕНИЕ УРАВНЕНИЙ КИНЕТИКИ ДВОЙНОЙ 
СВЯЗАННОЙ РЕАКТОРНОЙ СИСТЕМЫ ДЛЯ 
СКАЧКА РЕАКТИВНОСТИ 


Точное решение в аналитическом виде уравнений, 
установленных в предыдущей работе [1], можно получить 
только для одного случая-скачка реактивности. Здесь 
мы приведем и проанализируем решение для случая, ког- 
да реактивность вводится только в одну сборку. 


В дальнейшем нам понадобится условие критично- 
сти. Его легко получить из уравнений кинетики, если по- 
ЛОЖИТЬ 


6; — бо1, б> — Бьз, 
где в; И 25. — значения подкритичностей сборок, обес- 
печивающие точно критическое состояние системы. 


В этом случае уравнения кинетики сводятся к 
Вот Пе: — —- 112 Поз, (1) 
Роз Поз == -- ]21 Пл, (2) 


Где По, И П,› стационарные средние плотности нейтро- 
нов в сборках. Перемножая (1) на (2), получим условие 
критичности двойной связанной системы: 


Ро: Ро2 — Ти» Ти. (3) 


Пусть теперь в критическую систему вводится скачком 
положительная реактивностьр в области первой сборки, 
т. е. в уравнениях кинетики надо положить 


бт -= Ре1 Е Б, 
(4) 


2» — ВБь>. 


Начальные условия задаются соотношениями: . 


В. 3. 
| п, (+) И 
р № 


: 
В. 1 
п. (1) = По ДЛЯ < ты. 
Ё Г : 
Решение уравнений кинетики для условий (4) и (5) 
может быть получено в виде: 


| - Г. 
поль У’ А(рдехр@й,, 
} 





м Си (= п» Для 1< 0, 


(5) 








7 Ср (И Ее 


(6) 


> 


И -- 
п. ({) = п. > А. 191 ехр! 6), 
| 


где || — число групп запаздывающих нейтронов в пер- 
ВОЙ сборке, 
|› — то же во второй. В дальнейшем мы будем считать 
Г Ая [, — | 





Р) — есть простые корни характеристического урав-., 
нения: 
р ВН 
В — — В, ыы РЕНИ т 
.Р 1 = р -- ^; 
(7) 
и | рВ.> 
-- Ро Р0зехр [- (< ты | нь >. 
от 202 ехр [ — (112 - <> р] | эр — ТЕ, р | 


Заметим, что в (7) с помощью условия критичности 
(3) удалось избавиться от коэффициентов связи, так что 
в решениях `они не появятся. 


Далее, 


3. 
А, (р;) — и хе р. не —В, р,ехр [— (512 ад 
(8) - 


ВЕ: 
+5 р А Рот роз ехр [-- < г р; с, 
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где ы 


` 
. 


В | г у 8:2 | | '  руВь | 
== 2 Г в 2 р; — Вь> .. Е 
| = р; м | | ы | р; + ^, 


| . в. 
в р г 12 
С, == |»: ры ы 15 р;- а | ОЭ -— р + 








гри 
Г В р;В» \“ 
+(ь- У ,)} (ру - хе ей а 
рН арм 7 ° 
А | В; в | 
в) |5 - 2902 )  — 35 - 
2(р;) = | И р ати р роз | ехр! — тру 
(9) 
Г. АВ 
| т РЗ + м -- Ролро> ехр [-— (112 -Ё т.) РЯ С, | 
где 





Ва р; В 
В. = ПУ Пр; — в — В = |, 
2 | 1 Ра Ей | "Ру — бо — В р А 


С, = | (вать) (| — р. } + 
== 12 1 721 \Р: 1 — Р 7” р 


Ви | Г рвы \ 
| В Е 
+ ( + а Я} [в о-ь-р+ $, ие 


Если р›>0, т. е. система надкритична, то один из кор- 





1 
ней уравнения (7) Ро.= — будет положительным. Другие 


корни отрицательные и соответствующие им члены в ре- 


шении (6) будут затухать со временем. Всего корней 
2(1+1). 


Если р <0, то все корни отрицательные, наконец, 
если р=0, то один корень равев нулю, а все остальные 
отрицательные. 
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Сравнивая выражения (8) и (9) с аналогичным вы- 
ражением для коэффициентов решения при скачке реак- 


тивности в однозонном реакторе [2]: 
э 


ое | в в а 
В 
| еее РЕ А \ Е“ (Ру м)“/’ 
(10) 


легко заметить, что влияние другой сборки учитывается 
добавлением третьего члена в числителе и знаменателе 
формул (8) и (9). 

В решении. можно различить несколько стадий: 
«мгновенный скачок», переходный ‘режим и установив- 
шееся экспоненциальное изменение. Мгновенный скачок 
связан с наибольшим отрицательным корнем Р_ уравне- 
ния (7). | 

Могут встретиться две ситуации. 

(1). Если |Р |», что всегда выполняется, но |Р_ |< 


] { . 





О 
$135 - 
1 и бы - В — В, | 
ое Ее 
Р то роз — В. 
о _. бо? _ 1 — В, \ (1) 
ре бо> — 2 т те 


Для однозонного реактора аналогичное выражение 
имеет вид 


1 | 


— 


Р СРВ. 
] | 
2). Если |Р/_ . и |Р_ —, то пПро- 
(2). Если |Р-1» м и Р-|% р 
стого выражения для Риз трансцендентного уравнения 


(7) получить нельзя. 


После скачкообразного изменения реактивности от 0 
до р нейтронный поток за время —1/|Р_| изменяется 


ОТ По! ДО П, = Пу [1 гае А, (Р_)] И ОТ 
По? ДО П› = Поз [1 — А. (Р_)]. 
41 


После скачка наступает переходная стадия, при кото- 
рой отношения С; к п изменяются от величин, определяе- 
мых начальными условиями (5) и характеризующими 
стационарный режим при нулевой реактивности до вели- 
чин, соответствующих установившемуся режиму при ре- 
активности р. Можно считать, что переходные процессы 
затухают за время порядка нескольких установившихся 
периодов 1/Р.. 


Установившийся период оценивается формулой 
1 Со РЕ: по НЙ: | 
— м | — И -- 21 -- НЕ о 
‚Ви? | Во2 
(12). 
Для однозонного реактора аналогичное зыражение 
имеет вид В 





Здесь ив (12); В+{-- УЗ, ^...Т. о.. первая сборка при 

1 
вводе реактивности о ‘ведет себя так, как будто время 
жизни мгновенных нейтронов в этой сборке увеличилось 


на 1. В ‚ а время жизни запаздывающих нейтро- 


Во2 
ры В 


НОВ В ЭТОЙ сборке унеличилось на Ре я 
09 


Заметим еще, что стадия скачка и стадия переходно- 
го режима усложняются явлениями запаздывания. 

Так из (6) и (8) видно, что решение для первой сбор- 
ки состоит из двух частей: изменение потока, связанное 
со свойствами первой сборки, проявляется гразу же в 
момент возмущения и описывается первой суммой чле- 
нов правой части (8). 

Изменение, связанное с влиянием второй сборки, про- 
является через время (112-121). 

Поток во второй сборке, после «мгновенного» ввода 
реактивности в первую, изменяется с запозданием то! в 
сравнении с основным изменением в первой. 

Посколько в решениях для п!({) и п2({) в показате- 
лях экспонент стоят корни одного и того же характери- 


_—— 
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стического уравнения, то при 1» (т;2+т2:) изменения по- 
токов в обеих сборках будут происходить одинаковым 
образом. 
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Е. Д. ВОРОБЬЕВ, Е. Г. КАЛАШНИКОВ 


ОБЩЕЕ РЕШЕНИЕ УРАВНЕНИЙ КИНЕТИКИ 
ДВОЙНОЙ СВЯЗАННОЙ РЕАКТОРНОЙ СИСТЕМЫ 


Подавляющее большинство задач реакторной кинети- 
ки может быть решено только приближенными метода- 
ми. Здесь мы попытаемся найти решение уравнений ки- 
нетики двойной связанной системы [1] методом, который 
применялся Аккези [2] при анализе кинетики однозон- 
ного реактора. ` 

Для простоты ограничимся одной группой запазды- 
вающих нейтронов. Члены с запаздыванием разложим в 
степенной ряд 


п (Е — =) = п (4) — зп (+) -... 
и оборвем его на втором члене. Это, очевидно, можно сде- 
лать для {< т и даже для {= т, если 
п (1) «п(Ф. 


Тогда исходные уравнения будут иметь вид: 








р, —В 112 12 
по) = и, (9) + р п» (9 — торс, 
я аа р» — В. 12 21 ° 
п›(#) = - , п О-о — ыы (И АС. ©, 
(1) 
С, (= —^С, (+ ты п, (1), 
ИИ 





С. (9 = —^С, (9 + п» (1). 


1, 
Здесь точки означают дифференцирование по {. Пусть 
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в одной из сборок двойной системы, скажем первой, ре- 
активность меняется по произвольному закону: 


р1 = Ри, -| Ро -| 22 (1), а. 


ыы (2) 
р2 — Рь2. 
Здесь р» — постоянная часть избыточной реактив- 
ности в первой сборке; 
2 (1) — функция, характеризующая временную зави- 


симость реактивности, а 
р— амплитуда изменений реактивности, причем . 
дем пытаться искать решение в виде: 


и — ехр [то (+) 91 


{ 
С, (+) =Р, (бехр [т (+) а: ). 


{ 3 
пы) = (дер (11 (94 ), (9) 


{ 
С; (1) == Р, ( ехр | ш (7) 91 
О 
Подставим выражения. (3) в уравнение (1): 
м (9) = Ро! - Ро + 22(+) — В, -- 1124 (© — 2129) -- 
— 512 112 Ч (+) т (1) + Ар, (1), 
154 (4) + 124 (т(0 -- (оз а 8.) а(+) -- 121 — т 12: 1 (1) -- 
--А1Ь р (1), (4) 
1.р: (#) Е | ри) м (9 = — М,р: © + 3, 
1 рг (#) Е р (1) п (9 = — М р.) + В» 9(#) 
Положим, что т (4), 4({), Р! (+) и Р2({) можно разло- 
жить в ряд по степеням р: 
п1(+) = по -- ри: (Е рт) .. 
4({) = 4 р 4; (& + а. (®-..., | 
ра(+) = рло + р р:1 (4) + р,» 4) +..., (5) 
рэ(+) = рго + рр», (+) + р р» (®-.... 


Здесь. шт; (1), а; (1), Ри(#) и Рь ({) — ограниченные функ- 
ции времени. Это требование обеспечивает сходимость 
рядов. 

° Шо Рь, Р»ь и 9, — постоянные. 

Подставим разложения (5) в уравнения системы (4) 
и уравняем коэффициенты при одинаковых степенях р. 
Тогда для нулевой степени о имеем: 


То = би Е рь — В: | 11? 4 — Т12 112 Чо Шо + АН Рио, 
[5 = (Риз — В») Чо -- 7121 — ть Ти Шо -- ^ 15 Рь, 
по 1: Ри —= —^ В Рь В, (6) 
По! Ро = — А Рьь - В Чо. 


Откуда 
Я | В 1 
и а = 7 
Ро Г ША" Ру = 4% ‚А (7) 
а 
в | я ^ 3 ; (1 ) 
Че = | Ши: — Ви — 90 Е В, = 9 | — 312 Шо 


ИЛИ 


/ 


В. 
Чо = 1»! (1 — вить | ть — Ро т — ет) 


(8) 
Исключая 4о из (8)`и (8!), получим: 
го ` А (Шо) = 1, (9) 
где % 
ы __ В То 
А (ть) = той — ри + =. 


Шо 
‚-— [бо Ро> (1 > Шо) (1 АА п.) ] | По |. —_ Роз ш и 


„ (10) 
Уравнение (9) есть модификация уравнения обрат- 
ных часов для двойной системы [3], а (10) есть переда- 


точная функция для такой системы, как мы ее опреде- 
ляем в линейной теории. 
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Для первой степени имеем: 


На (И) = 214) 112, (#) — 1291 (И) — 11а (Й — 
— р 2 Чи Ш: в НА Ри, (+), 


1, Ч (Им, 9, (©) — 1. Че т; (1) = (5,5 — В} 4: (#} — 


— та Ти Пи (И -= 2.1 р», (1, (11) 
ри (#) им, ри (4) -- п рьм, (6 = — М. ри (4), 


Тр (4) + тьрь: (1) -- Ро п; (#} = — М.рь(® - 8,99. 


Эти уравнения можно проинтегрировать с помощью 
преобразований Лапласа. Членами, содержащими т Т, 
будем в дальнейшем пренебрегать. Можно показать, что 
они меньше других ‘членов уравнений (11). 


им, ($) == 2($) + 11-9 $) -Е М:ри (3), 
1.391(3) —- 15 ть 41(3) - 15 Ч. м3) = 
— (602 — Вэ) 4: ($) + А 1. р», ($), (12) 
1 Зри, (3) ++ Ито ри: ($) -- приют, ($) = — М р1:($), 

15$ рэ! ($) - 1ьто ры ($) Е 15рэь 11 ($) = — Мэри ($) 8.43). 

Здесь п, ($), р11($),рь ($) и 4.($) изображения функций 
т; (+), риг (+), р»! (+) и Чт (1), соответственно. 

После несложных преобразований получим: 


п: ($) == К($ па,) . 2 (3), (13). 
где 


{= ыы о а а 
а ее 


Ро в Г, } Шь -- ^ $ Не По ть х)/ 
и 1. 
| 15 ($ Ее а) АР 35 = 5 -- п. -- ] 
| По 3 у 
О аитий ао ——- ыы ы 14 
| баз = по + | (14) 
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Заметим, что 





АВ 
м а: а аа 
К (0: По) ве е! не (110 а ). 2 ря 
| ^В, Шов» \” 
— Рог Виз [в + Риты + га | 
(14') 
а К-Цо: о) = 1.1 ы (141) 
02 
где 
Бен, 
а В = + _- 


Подобным же путем можно получить выражения для 
12 (1) ит. д., но вычисления, к сожалению, становятся 
очень громоздкими и рассматриваться здесь не будут. 

Рассмотрим теперь несколько конкретных примеров. 


1. Ступенчатая функция 


Если в уравнении (2) р=О0, а р-=0О — получаем 
скачок реактивности и с помощью уравнения (9) полу- 
чаем точное решение для этого случая: 


п; (9 = № С; ехр(ть®), (15) 


как в работе [3]. Решение можно ‘представить и в 
другом виде. Пусть ро =0, 2(1) =1. Можно показать [4], 
что 

т, (#) =К (01 11%}. 


Тогда, в первом приближении, общее решение будет 
п; (В) = С. ехр (рК (0,0)1+) (16) 
-- > Сехр[ты Е РК (01 По!) |1}, (17) 


где Сои С; — постоянные интегрирования, которые дол- 
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жны быть определены из начальных условий для 
п; (Е), п2({) ‚С, (® иС2(), ам, — отрицательные корни 
уравнения (9). 

Можно заметить, что определяющим в решении (17) 
будет первый член. Т. о., в установившемся режиме первая 
сборка ведет себя так, как будто общее время жизни 
нейтронов, согласно (14!) увеличилось на В. 

03 


2. Линейная функция 


Если предыдущую задачу можно точно решить пря- 
мым методом и ее решение приведено здесь лишь для 
сравнения, то эта задача не тривиальна. 


Пусть 


Визбыт. — Во нЕ РЁ, 
т. е. 2 (1) =. 


В первом приближении решение может быть получе- 


но из (13) подстановкой 2 (3) = 1/3?: 


(0 =УС, ехР] мые -- СРК (©: по) 8 | | (18) 


2 


здесь т‚; — все (четыре) корни уравнения (9). 

Очевидно, это уравнение справедливо в течение огра- 
ниченного промежутка времени, поскольку 2({) не огра- 
ничено во времени. 


3. Синусоидальная функция 


ехр (]е4) 
о 


Физический смысл имеет только действительная 
часть этого выражения. 


Пусть Ризбыт: — Во + р 


Тогда т,($) =— К ($115) к Кб 


Отсюда 
т, (+) = —К (4, По) | ехр (] 61) 
49 


| | 
и п, (0 = ГК (ель) рехр (494) @&= 
О (©) 


- К (9, шо\ехр (11). 


Введем обозначение 





® [< (1%, т.) -- 2 (]®, По). 


Тогда 
пи = © Сехр иьй 9 [2 6, пе) | $1 (в + +) |, 
где $; фазовый угол 2 (]®, Шо). 


Если система критична, т =0, а ш, дают зату- 
хающие члены и 


п; (1) = С,ехр [в |2(]®) | 51п («Е -[ $о)]. (19) 
При р/2(])| < | | 

п: (0 = С [1 +22 (}=)| Зи («Е - $о)]. (19') 
Здесь | 2, (]®)| играет роль передаточной функции. 
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